High momentum response of liquid 3 He 
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A final-state-effects formalism suitable to analyze the high-momentum response of Fermi liquids is 
presented and used to study the dynamic structure function of liquid 3 He. The theory, developed as 
a natural extension of the Gersch-Rodriguez formalism, incorporates the Fermi statistics explicitely 
through a new additive term which depends on the semi-diagonal two-body density matrix. The 
use of a realistic momentum distribution, calculated using the diffusion Monte Carlo method, and 
the inclusion of this additive correction allows for a good agreement with available deep-inelastic 
neutron scattering data. 
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Inelastic neutron scattering is the most efficient tool to 
explore the structure and dynamics of quantum liquids 
4 He and 3 He since the dynamic structure function S(q, to) 
is readily obtained from the double differential scattering 
cross-section Q . The range of momenta q transferred to 
the system determines the kind of microscopic informa- 
tion that can be extracted. The most interesting regimes 
correspond to low and high q's. At low q, the scatter- 
ing data allows for the determination of the low-energy 
excitation spectrum. In the opposite limit, known as 
deep-inelastic neutron scattering (DINS), q is so high 
that single-particle properties of the system become ac- 
cessible. 

It is well known that in the q — > oo limit, S(q,uj) ap- 
proaches the impulse approximation (IA). The only in- 
gredient to calculate the response in IA is the momentum 
distribution n(k), a fundamental function in the study 
of 4 He, 3 He, and the 4 He- 3 He mixture. The boson and 
fermion quantum statistics of 4 He and 3 He, respectively, 
introduce significant differences in their corresponding 
momentum distributions. Liquid 4 He presents a macro- 
scopic occupation of the zero-momentum state, charac- 
terized by its condensate fraction no; n{k) of liquid 3 Hc, 
considered as a normal Fermi liquid, shows a discontinu- 
ity at the Fermi momentum kp Q- Nowadays, different 
theoretical calculations of n(k) ranging from variational 
theory, based on the (Fermi-)hypernetted-chain equa- 
tions ((F)HNC), to the more exact diffusion Monte Carlo 
(DMC) method are in an overall quantitative agreement 
HiQild- However, a direct comparison with experimen- 
tal data is not possible due to instrumental resolution 
effects (IRE) and, more fundamentally, to final state ef- 
fects (FSE). From the theoretical side, the problem is 
that the IA does not account completely for the scat- 
tering in most of the DINS experiments since the trans- 



ferred momenta are not high enough. Therefore, FSE 
which take into account the interactions of the struck 
atom with the medium can not be disregarded. 

The search for an unambiguous experimental signa- 
ture of no in liquid 4 He using DINS has originated a 
great deal of theoretical and experimental work for the 
last two decades. At present, theoretical predictions 

H HSU [13 for both thc FSE and thc IA P rovide a sat " 

isfactory description of the experimental measurements, 
with an overall agreement on the value of the condensate 
fraction, no ~ 9% at the equilibrium density. It is worth 
noticing that FSE in superfluid 4 He are enhanced due to 
no and therefore, even at the highest momenta achieved 
in the laboratory, FSE play a fundamental role. Com- 
paratively, few works have been devoted to the analysis 
of the high-g response of liquid 3 He. The main reasons 
underlying this situation are, from the experimental side, 
the large neutron absorption of 3 He, and from the the- 
oretical one, the difficulties the Fermi statistics of 3 He 
introduces in the quantum many-body calculations. The 
most accurate data have been reported by Azuah et al. 

, and more recently by Senesi et al. 0| > but only the 
first one was carried out at the equilibrium density. FSJ3 
in 3 He have been taken into account by Moroni et al. 
using the bosonic formalism of Carraro and Koonin 
Their results || show less strength at the peak than the 
experimental S(q,u)), pointing to possible limitations of 
the formalism when applied to a Fermi liquid. On the 
other hand, an analysis of the experimental data based 
on cumulant expansions 0, ^| has revealed significant 
differences between the experimental and theoretical He 
momentum distributions at equilibrium density 

We present in this letter results for the high-g re- 
sponse of 3 He using a theoretical formalism that incor- 
porates explicitly and consistently the Fermi statistics to 
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the FSE. The inputs required are the momentum dis- 
tribution n(k) and the semi-diagonal two-body density 
matrix p 2 {r\, r 2 ; r[, r 2 ). Both are obtained from micro- 
scopic theory: n(k) from a DMC calculation, and p 2 from 
variational FHNC. The results obtained for S(q, to) repro- 
duce the experimental data better than previous estima- 
tions, pointing to non-negligible Fermi contributions to 
the FSE. 

As long as 4 He is concerned, most theories introduce 
FSE as a convolution in energies, which turn into an al- 
gebraic product in time representation. Hence, FSE are 
included by means of a new function R{q, t) which mul- 
tiplies Sia(<7, i) to obtain the total response, S(q,t) = 
Sia(<?, t)R(q, i). This convolutive approach is clearly am- 
biguous when applied to fermions because, contrary to 
the Bose case, >Sta(q, t) has an infinite number of nodes, 
a fact that leads to a singular definition of R(q,t). In 
order to overcome this serious drawback that emerges 
from a direct translation of the FSE theories for bosons 
to fermions, we have used an alternative formulation 
that can be considered a natural extension of Gersch- 
Rodriguez theory [LH to fermionic systems. In previous 
works 0, ^| , we have applied this theory to evaluate 
FSE in 4 He- 3 He mixtures, but there the fermionic cor- 
rections are much smaller due to the low 3 He concentra- 
tions. 

At high q, the density-density correlation factor S(q, t) 
can be well approximated by 

I j '^Q_ s f 
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with io q — q 2 /2m, v q — q/m, s = v q t, and AV(r,u) = 
V(r — u) — V(r). Equation is derived assuming that 
the atom struck by the neutron recoils in a medium of 
non-moving 3 He atoms. S(q,t) is still hard to evaluate 
using Eq. (yi since it implies an integration over the com- 
plete semi-diagonal ./V-body density matrix of the sys- 
tem, which is essentially the square of the ground-state 
wave function. In He, a truncated cumulant expansion 
of Eq. leads to the Gersch-Rodriguez expression for 
the FSE [l4|. Recently, this formalism has proven its 
efficiency by reproducing 4 He DINS data with high ac- 
curacy 7]. On the contrary, the nodal structure of pn in 
a Fermi system prevents from a straightforward applica- 
tion of these methods. 

In order to extend the FSE theory to 3 He, one in- 
troduces an auxiliary ./V-body density matrix p^ corre- 
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sponding to a system of spinless bosons with the mass, 
density, and interatomic potential of 3 He. p^ is positive 
and then it can be used as the starting point of a cumu- 
lant expansion of the response function. In terms of pfj, 
a convenient decomposition of pjy turns out to be 



p N (rx, . . -,r N ;r[) = pi(r n >) 



1 



pfVii') 



xpf I (r 1 , . . .,r N ;r[) 



+ Ap N (r 1 ,...,r N ;r[) , (2) 



Pi{fii') being the one-body density matrix extracted 
from p^. In the thermodynamic limit, pf (riv) factorizes 
from p^r, and thus the first term in Eq. (J2J) corresponds to 
an artificial TV-body density matrix containing fermionic 
correlations between points 1 and 1' only. 

Inserting p^ © in S(q,t) the 3 He response be- 
comes 



S(q,t) = S IA (q,t)R(q, t) + AS(q,t) , 



(3) 



with SiA(q,t) the exact 3 He IA, and R(q,t) the Gersch- 
Rodriguez FSE function calculated with the bosonic 
semi-diagonal two-body density matrix (t*i, r 2 ; r[), 



R(q, t) — exp 
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The new additive term AS(q, t) in Eq. © is a con- 
sequence of Apjv introduced in Eq. J2J). The leading 
contribution to the FSE at high q depends on the semi- 
diagonal two-body density matrix, and then Ap 2 is re- 
quired for the calculation of AS(q,t). A cluster expan- 
sion in the framework of the FHNC formalism allows for 
an estimation of Ap 2 according to the following struc- 
ture, 

-Ap 2 (ri,r 2 ;ri) = pi{r xv )G{r u r 2 ; r[) 
P 

+ pMnv)F(r 1 ,r 2 ;r' 1 ) . (5) 

The form factors G(ri,r 2 ;r[) and F(ri,r 2 ;r[) can be 
expressed in terms of auxiliary functions defined in the 
FHNC theory, and /Oid(^ii') is positive everywhere and 
similar to a bosonic one-body density matrix [171 Il8j . 
Therefore, Pid(^h') can be used as the basis of a cumu- 
lant expansion by simply adding and subtracting it to 
ApN. The resulting additive correction AS(q, t) is, to 
the lowest order, 



AS(q,t)=e iu " s ^-p w (r lv ) 



exp 



3_ f drAp 2 (r,0;r + s) 

PiT>{s) J 



1 — exp 



ds'AV{r,s'] 



(6) 
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Finally, a Fourier transform of S(q, t) provides the dy- 
namic structure function S(q,u>). Furthermore, the scal- 
ing properties of the IA, in terms of the West variable 0] 
Y = mui/q — q/2, suggests as usual to write the response 
{q/m)S{q 1 u)) as a Compton profile, 

J(q, Y) = J dY'J(Y')R(q, Y-Y') + AJ(q, Y) , (7) 

J(Y) = 1/(2tt 2 p) C, dkkn(k) being the IA. 

The microscopic functions entering the high- 
momentum response J(q, Y) (JZJ are the one- and 
the semi-diagonal two-body density matrices of the 
actual system and its bosonic counterparts. The most 
relevant quantity is the one-body density matrix, or 
equivalently the momentum distribution, which is used 
to evaluate J(Y). We have estimated the 3 He momen- 
tum distribution using the DMC methodology that has 
recently proved to be very accurate in the calculation of 
the 3 He equation of state at zero temperature |2fJ, |2l| . 
At the equilibrium density po = 0.273 ct~ 3 (<t = 2.556 
A), n{k) is well parameterized by 

( a Q - a 3 x 3 x < 1 

n[X> ~\{bo + hx + b 2 x 2 )e- b ' x x > 1 ' [ ) 

with x = k/kp, k-p being the Fermi momentum. The set 
of parameters that best fit n(k) (JHJ) is reported in Table 
I. The kinetic energy per particle, related to the second 
moment of n(k), is 12.3 K and the discontinuity of n{k) 
at the Fermi surface is Z = 0.236. The value of Z, which 
defines the strength of the quasi-particle pole, is rather 
small, indicating that the system is strongly correlated. 
On the other hand, the tail of the momentum distribution 
extends up to high momenta generating significant high- 
energy wings in J(q, Y). The present n(k) is in overall 
agreement with the DMC one from Ref. [3j. 

The semi-diagonal two-body density matrix p 2 , and 
the auxiliary bosonic functions pf and p^ , have been ob- 
tained in the framework of the FHNC and HNC theories 
using a Jastrow-Slater variational wave function. It is 
well known that this trial wave function is not accurate 
enough if the main objective is to get a good upper-bound 
to the energy. This is not however the aim of the present 
letter. In fact, we have shown in previous work that a Jas- 
trow wave function can efficiently account for the FSE in 



a 0.481319 

a 3 0.0842956 

b 1.39056 

bi 0.157930 

b 2 0.0829832 

b t 2.31398 



TABLE I: Parameters of n(k) JHJ at po- 



4 He 0. Certainly, short-range correlations, which domi- 
nate the FSE, are already contained in the Jastrow-Slater 
approximation. Accordingly, the diagrammatic analysis 
for 3 He has been performed at the two-body level, thus 
making the analysis easier. 

The FSE broadening function R(q, Y) at p and a mo- 
mentum transfer q — 19.4 A -1 is shown in Fig. ^ This 
value of q has been used throughout this work since it cor- 
responds to the momentum reported in the experimental 
data by Azuah et al. ^(?>^) has been calculated 

in the bosonic approximation and then its structure is 
similar to the FSE function of 4 He [7j. When q increases 
R{q, Y) narrows and sharpens, becoming a delta function 
in the q — > oo limit. 

The 3 He additive correction A J(q, Y) = 
(q/m)AS(q, Y) at the same density and momen- 
tum transfer is also shown in Fig. ^ This function, 
which introduces fcrmionic correlations to the FSE, 
presents a shape that is entirely different from that 
of R(q,Y). The strong oscillations that appear in the 
region Y w ±fcp modify the shape of the IA response 
around these points. Furthermore, a central peak 
centered at Y = enhances the strength of the total 
response at the origin. Out of this region (\Y\ > /cf), 
A J(q, Y) is much smaller and its correction to the 
response becomes negligible. Further analysis indicates 
that A J(q, Y) decays to zero in the high momentum 
transfer limit. This fact, together with the limiting 
condition R{q — > oo,y) — > S(Y), indicates that the total 
response asymptotically approaches J(Y) when q — > oo. 

The final result of the 3 He response at q = 19.4 A -1 
is shown in Fig. [3 We compare our results with the 
DINS data of Azuah et al. ,llj because their data cor- 
respond to densities around the equilibrium density po, 
and also because their experimental setup produces a 
rather narrow instrumental resolution function. More re- 
cent data 0| are focused at higher liquid densities and 
to the solid phase. In this experiment, the momentum 
transfer is much larger (q ~ 90 to 120 A -1 ) produc- 
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FIG. 1: FSE broadening function R(q, Y) (solid line) and 
additive correction AJ(q,Y) (dashed line) at q = 19.4 A -1 
and po- AJ(q, Y) is multiplied by a factor 20 to fit into the 
scale. 
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FIG. 2: 3 He response at q=19.4 A -1 and po (solid line), folded 
to the IRE function (dotted line). Thepoints with error bars 
are the experimental data from Ref. pl|l . The instrumental 
resolution function [TH has been divided by a factor 10 to fit 
into the scale. The inset shows J(q, Y) near the peak with 
(solid line) and without (dashed line) the additive correction 
AJ(q,Y). 

ing a simultaneous decrease of the FSE corrections and 
a widening of the IRE. The IRE function estimated by 
Azuah et al. ^lj is shown in Fig. [21 scaled by a factor 
0.1. The theoretical response J(q, Y), shown in the fig- 
ure, has been folded with the experimental IRE function 
I(q, Y) to make a direct comparison possible. The kinks 
of J(q, Y) at Y = ±kp present in the IA, are completely 
washed out by the succesive folding with R(q,Y) and 
I(q,Y), although AJ(q,Y) introduces additional struc- 
ture around those points. The most remarkable feature is 
the enhancement of the strength of the response around 
the peak due to the new addititve term introduced in the 
present approach (see inset in Fig. [2J. This small but 
significant increase of strength allows for the first time to 
reproduce the available 3 He DINS experimental data. 

To summarize, we have presented a FSE formalism 
suitable to study the dynamic structure function of a 
Fermi liquid like liquid 3 He at large momentum trans- 
fer. The method is a natural extension of the Gersch- 
Rodriguez theory for bosons. According to the present 
formulation, J(q, Y) results from the convolution of the 
IA with a purely bosonic FSE broadening function, which 
incorporates the short-range correlations induced by the 
interatomic potential, plus an additive correction term 
that takes into account Fermi statistics effects in the FSE. 
The results obtained are in good agreement with DINS 
data, comparable for the first time to the accuracy pre- 



viously achieved in the study of the high-q response of 
liquid 4 He. The two key features behind the present re- 
sults are, on one hand, the use of a realistic momentum 
distribution provided by the DMC method, and on the 
other, the explicit introduction of Fermi corrections in 
the FSE. The latter effect is estimated at the lowest or- 
der but its inclusion allows for a significant improvement 
and a better knowledge of specific mechanisms influenc- 
ing the FSE in liquid 3 He. 
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